A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements

Study of Aerothermodynamic Modeling Issues Relevant to High-Speed Sample Return Vehicles

This paper examines the application of state-of-the-art coupled ablation and radiation simulations to highspeed sample return vehicles, such as those returning from Mars or an asteroid. A defining characteristic of these entries is that the surface recession rates and temperatures are driven by nonequilibrium convective and radiative heating through a boundary layer with significant surface blowing and ablation products. Measurements relevant to validating the simulation of these phenomena are reviewed and the Stardust entry is identified as providing the best relevant measurements. A coupled ablation and radiation flowfield analysis is presented that implements a finite-rate surface chemistry model. Comparisons between this finite-rate model and a equilibrium ablation model show that, while good agreement is seen for diffusion-limited oxidation cases, the finite-rate model predicts up to 50% lower char rates than the equilibrium model at sublimation conditions. Both the equilibrium and finite rate models predict significant negative mass flux at the surface due to sublimation of atomic carbon. A sensitivity analysis to flowfield and surface chemistry rates show that, for a sample return capsule at 10, 12, and 14 km/s, the sublimation rates for C and C3 provide the largest changes to the convective flux, radiative flux, and char rate. A parametric uncertainty analysis of the radiative heating due to radiation modeling parameters indicates uncertainties ranging from 27% at 10 km/s to 36% at 14 km/s. Applying the developed coupled analysis to the Stardust entry results in temperatures within 10% of those inferred from observations, and final recession values within 20% of measurements, which improves upon the 60% over-prediction at the stagnation point obtained through an uncoupled analysis. Emission from CN Violet is shown to be over-predicted by nearly and order-of-magnitude, which is consistent with the results of previous independent analyses. Finally, the coupled analysis is applied to a 14 km/s Earth entry representative of a Mars sample return. Although the radiative heating provides a larger fraction of the total heating, the influence of ablation and radiation on the flowfield are shown to be similar to Stardust

Uncertainty Analysis of Air Radiation for Lunar Return Shock Layers

By leveraging a new uncertainty markup technique, two risk analysis methods are used to compute the uncertainty of lunar-return shock layer radiation predicted by the High temperature Aerothermodynamic Radiation Algorithm (HARA). The effects of epistemic uncertainty, or uncertainty due to a lack of knowledge, is considered for the following modeling parameters: atomic line oscillator strengths, atomic line Stark broadening widths, atomic photoionization cross sections, negative ion photodetachment cross sections, molecular bands oscillator strengths, and electron impact excitation rates. First, a simplified shock layer problem consisting of two constant-property equilibrium layers is considered. The results of this simplified problem show that the atomic nitrogen oscillator strengths and Stark broadening widths in both the vacuum ultraviolet and infrared spectral regions, along with the negative ion continuum, are the dominant uncertainty contributors. Next, three variable property stagnation-line shock layer cases are analyzed: a typical lunar return case and two Fire II cases. For the near-equilibrium lunar return and Fire 1643-second cases, the resulting uncertainties are very similar to the simplified case. Conversely, the relatively nonequilibrium 1636-second case shows significantly larger influence from electron impact excitation rates of both atoms and molecules. For all cases, the total uncertainty in radiative heat flux to the wall due to epistemic uncertainty in modeling parameters is 30% as opposed to the erroneously-small uncertainty levels (plus or minus 6%) found when treating model parameter uncertainties as aleatory (due to chance) instead of epistemic (due to lack of knowledge)

A Boundary Condition Relaxation Algorithm for Strongly Coupled, Ablating Flows Including Shape Change

Implementations of a model for equilibrium, steady-state ablation boundary conditions are tested for the purpose of providing strong coupling with a hypersonic flow solver. The objective is to remove correction factors or film cooling approximations that are usually applied in coupled implementations of the flow solver and the ablation response. Three test cases are considered - the IRV-2, the Galileo probe, and a notional slender, blunted cone launched at 10 km/s from the Earth's surface. A successive substitution is employed and the order of succession is varied as a function of surface temperature to obtain converged solutions. The implementation is tested on a specified trajectory for the IRV-2 to compute shape change under the approximation of steady-state ablation. Issues associated with stability of the shape change algorithm caused by explicit time step limits are also discussed

Features of Afterbody Radiative Heating for Earth Entry

Radiative heating is identified as a major contributor to afterbody heating for Earth entry capsules at velocities above 10 km/s. Because of rate-limited electron-ion recombination processes, a large fraction of the electronically-excited N and O atoms produced in the high temperature/pressure forebody remain as they expand into the afterbody region, which results in significant afterbody radiation. Large radiative heating sensitivities to electron-impact ionization rates and escape factors are identified. Ablation products from a forebody ablator are shown to increase the afterbody radiation by as much as 40%. The tangent-slab radiation transport approach is shown to over-predict the radiative flux by as much as 40% in the afterbody, therefore making the more computationally expensive ray-tracing approach necessary for accurate radiative flux predictions. For the Stardust entry, the afterbody radiation is predicted to be nearly twice as large as the convective heating during the peak heating phase of the trajectory. Comparisons between simulations and the Stardust Echelle observation measurements, which are shown to be dominated by afterbody emission, indicate agreement within 20% for various N and O lines. Similarly, calorimeter measurements from the Fire II experiment are identified as a source of validation data for afterbody radiation. For the afterbody calorimeter measurement closest to the forebody, which experiences the largest afterbody radiative heating component, the convective heating alone is shown to under-predict the measurement, even for the fullycatalytic assumption. Agreement with the measurements is improved with the addition of afterbody radiation. These comparisons with Stardust and Fire II measurements provide validation that the significant afterbody radiation values proposed in this work are legitimate

Impact of State-Specific Flowfield Modeling on Atomic Nitrogen Radiation

A hypersonic flowfield model that treats electronic levels of the dominant afterbody radiator, N, as individual species is presented. This model allows electron-ion recombination rate and two-temperature modeling improvements, the latter which are shown to decrease afterbody radiative heating by up to 30%. This increase is primarily due to the addition of the electron-impact-excitation energy-exchange term to the energy equation governing the vibrational-electronic-electron temperature. This model also allows the validity of the often applied quasi-steady state (QSS) approximation to be assessed. The QSS approximation is shown to fail throughout most of the afterbody region for lower electronic states, although this impacts the radiative intensity reaching the surface by less than 15%. By computing the electronic state populations of N within the flowfield solver, instead of through the QSS approximation in the radiation solver, the coupling of nonlocal radiative transition rates to the species continuity equations becomes feasible. Implementation of this higher- fidelity level of coupling between the flowfield and radiation solvers is shown to increase the afterbody radiation by up to 50% relative to the conventional model

Impact of Coupled Radiation and Ablation on the Aerothermodynamics of Meteor Entries

A high-fidelity approach for simulating the aerothermodynamic environments of meteor entries is developed. Two primary components of this model are coupled radiation and coupled ablation. Coupled radiation accounts for the impact of radiation on the flow field energy equations, while coupled ablation explicitly models the injection of ablation products within the flow field and radiation simulations. For a meteoroid with a velocity of 20 km/s, coupled radiation reduces the stagnation point radiative heating by over 60%. For altitudes below 40 km, the impact of coupled radiation on the flow field structure is shown to be fundamentally different, as a result of the large optical thicknesses, than that seen for reentry vehicles, which do not reach such altitudes at velocities greater than 10 km/s. The impact of coupled ablation (with coupled radiation) is shown to provide at least a 70% reduction in the radiative heating relative to the coupled-radiation-only cases. This large reduction is partially the result of the low ionization energies, relative to air species, of ablation products. The low ionization energies of ablation products, such as Mg and Ca, provide strong photoionization and atomic line absorption in regions of the spectrum that air species do not. MgO and CaO are also shown to provide significant absorption. Turbulence is shown to impact the distribution of ablation products through the shock- layer, which results in up to a 100% increase in the radiative heating downstream of the stagnation point. To create a database of heat transfer coefficients the developed model was applied to a range of cases. This database considered velocities ranging from 14 to 20 km/s, altitudes ranging from 20 to 50 km, and nose radii ranging from 1 to 100 m. The heat transfer coefficients from these simulations are below 0.045 for the range of cases (with turbulence), which is significantly lower than the canonical value of 0.1

A Study of Ablation-Flowfield Coupling Relevant to the Orion Heatshield

The coupled interaction between an ablating surface and the surrounding aerothermal environment is studied. An equilibrium ablation model is coupled to the LAURA flowfield solver, which allows the char ablation rate (m-dot(sub c)) to be computed as part of the flowfield solution. The wall temperature (T(sub w)) and pyrolysis ablation rate (m-dot(sub g)) may be specified by the user, obtained from the steady-state ablation approximation, or computed from a a material response code. A 32 species thermochemical nonequilibrium flowfield model is applied, which permits the treatment of C, H, O, N, and Si containing species. Coupled ablation cases relevant to the Orion heatshield are studied. These consist of diffusion-limited oxidation cases with Avcoat as the ablation material. The m-dot(sub c) values predicted from the developed coupled ablation analysis were compared with those obtained from a typical uncoupled ablation analysis. The coupled results were found to be as much as 50% greater than the uncoupled values. This is shown to be a result of the cumulative effect of the two fundamental approximations inherent in the uncoupled analysis

Nonequilibrium Stagnation-Line Radiative Heating for Fire II

This paper presents a detailed analysis of the shock-layer radiative heating to the Fire II vehicle using a new air radiation model and a viscous shock-layer flowfield model. This new air radiation model contains the most up-to-date properties for modeling the atomic-line, atomic photoionization, molecular band, and non-Boltzmann processes. The applied viscous shock-layer flowfield analysis contains the same thermophysical properties and nonequilibrium models as the LAURA Navier-Stokes code. Radiation-flowfield coupling, or radiation cooling, is accounted for in detail in this study. It is shown to reduce the radiative heating by about 30% for the peak radiative heating points, while reducing the convective heating only slightly. A detailed review of past Fire II radiative heating studies is presented. It is observed that the scatter in the radiation predicted by these past studies is mostly a result of the different flowfield chemistry models and the treatment of the electronic state populations. The present predictions provide, on average throughout the trajectory, a better comparison with Fire II flight data than any previous study. The magnitude of the vacuum ultraviolet (VUV) contribution to the radiative flux is estimated from the calorimeter measurements. This is achieved using the radiometer measurements and the predicted convective heating. The VUV radiation predicted by the present model agrees well with the VUV contribution inferred from the Fire II calorimeter measurement, although only when radiation-flowfield coupling is accounted for. This agreement provides evidence that the present model accurately models the VUV radiation, which is shown to contribute significantly to the Fire II radiative heating

Extension of Multiband Opacity-Binning to Molecular, Non-Boltzmann Shock Layer Radiation

For accurate predictions of shock layer radiative heating to reentry vehicles, the smeared rotational band (SRB) model is appropriate for molecular band systems with negligible self absorption, meaning they are optically-thin. However, for band systems with noticeable self absorption, the orders-of-magnitude more computationally expensive line-by-line (LBL) approach is required. Considering past and proposed NASA missions, the molecular band systems most likely to require the LBL approach are the CO 4th-Positive, CN Violet, and CO2 IR bands. The CO 4th- Positive and CN Violet bands are required for Mars entry at velocities greater than 6 km/s, with the CN Violet band also required for Titan entry. These two bands typically emit strongly in flow regimes with non-Boltzmann upper electronic state populations. The CO2 IR band is required for Mars entry at velocities below 5 km/s. This ro-vibrational band system is typically assumed to contain Boltzmann populations of radiating levels (the quality of this assumption is the subject of other studies)